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ABSTRACT 

Context. More than 60 % of Cepheids are in binary or multiple systems. Studying such systems could lead to a better understanding 
of the age and evolution of Cepheids. These are also useful tools to estimate the mass of Cepheids, and constrain theoretical models 
of their pulsation and evolution. 

Aims. We aim at determining the masses of Cepheids in binary systems, as well as their geometric distances and the flux contribution 
of the companions. The combination of interferometry with spectroscopy will offer a unique and independent estimate of the Cepheid 

masses. 

Methods. Using long-baseline interferometry at visible and infrared wavelengths, it is possible to spatially resolve binary systems 
containing a Cepheid down to milliarcsecond separations. Based on the resulting visual orbit and radial velocities, we can then derive 
the fundamental parameters of these systems, particularly the masses of the components and the geometric distance. We therefore 
performed interferometric observations of the first-overtone mode Cepheid VI 334 Cyg with the CHARA/MIRC combiner. 
Results. We report the first detection of a Cepheid companion using long-baseline interferometry. We detect the signature of a 
companion orbiting V1334 Cyg at two epochs. We measure a flux ratio between the companion and the Cepheid / = 3.10 ± 0.08 %, 
giving an apparent magnitude m H = 8.47 ± 0.15 mag. The combination of interferometric and spectroscopic data have enabled the 
unique determination of the orbital elements. P = 1938.6 ± 1.2 days, T p = 2 443 616.1 ± 7.3, a = 8.54 ± 0.51 mas, i = 124.7 ± 1.8°, 
e = 0.190 ± 0.013, a) = 228.7 ± 1.6°, and CI = 206.3 ± 9.4°. We derive a minimal distance d ~ 691 pc, a minimum mass for both 
stars of 3.6 M G , with a spectral type earlier than B5.5V for the companion star. Our measured flux ratio suggests that radial velocity 
detection of the companion using spectroscopy is within reach, and would provide an orbital parallax and model-free masses. 

Key words, techniques: interferometric - techniques: high angular resolution - stars: variables: Cepheids - star: binaries: close 



1. Introduction 

Classical Cepheid stars have been considered as reliable tools 
to estimate distances in the universe for more than a century 
(see e.g. Leavitt & Pickering 1912; Fernie 1969; Sandage & 
Tammann 2006; Barnes 2009; Bono et al. 2010). Their Period- 
Luminosity (P-L) relation makes them valuable to determine 
extragalactic distances and to calibrate secondary distance indi- 
cators. In addition to the determination of distances, Cepheids 
are also powerful astrophysical laboratories that provide fun- 
damental clues for studying the pulsation and evolution of 
intermediate-mass stars (see e.g. Prada Moroni et al. 2012; Bono 
et al. 2006; Caputo et al. 2005). 

The occurrence of Cepheids in binary (multiple) systems 
seems to be as high as 60 % for the brightest Cepheids (Szabados 
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2003), and is often neglected in the determination of the P-L re- 
lation. The position of a Cepheid in the P-L diagram could be 
biased by the presence of close bright companions. If the dif- 
ference in magnitude between the Cepheid and its companion is 
small, the apparent magnitude of the Cepheid will be overesti- 
mated. This has a particular impact on the use of a surface bright- 
ness (SB) method to estimate the radius and luminosity (e.g. 
Fouque & Gieren 1997; Gieren et al. 1998). The radial veloc- 
ity measurements can also be altered because of orbital effects. 
This leads to a bias in the distance estimate, since radial veloci- 
ties are necessary in the Baade-Wesselink method to evaluate the 
radius. The knowledge of potential photometric and astrometric 
biases caused by the companions is therefore important for the 
distance scale. 

Cepheids belonging to binary systems also offer the unique 
opportunity to make progress in resolving the Cepheid mass 
problem. For many years, stellar evolutionary models have pre- 
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Table 1. Parameters of the Cepheid V1334 Cyg and its close companion. 



Primary (Cepheid) 










Secondary^ 








m w " m K b m H b Sp. Type 1 ' P pul c 


Old" 


d e 


Sp. Type 


Porb 


To 


e ai sin i 


CO 


KM) 


(days) 


(mas) 


(pc) 




(days) 


(days) 


(AU) 


(rad) 


(Mo) 


5.87 4.46 4.66 F2Ib 3.333 


0.534 


683 




1937.5 


2 443 607 


0.197 2.46 


3.95 


0.529 



Notes. m v , m K , m H : mean apparent V, K and H magnitudes. Sp. Type: spectral type. P pul : period of pulsation. 9 LD : mean angular diameter, d: 
distance. P or t,: orbital period. 7b: time passage through periastron. e: eccentricity of the orbit. a\ sin;: projected semi-major axis of the orbit of the 
Cepheid about the center of mass of the system, to: argument of periastron. f(M): spectroscopic mass function. 

(fl) from Klagyivik & Szabados (2009). (A) from the 2MASS catalogue (Cutri et al. 2003). (c) from Samus et al. (2009). W) from the surface brightness 
relation of Kervella et al. (2004a) (f) from the P-L relation of Bono et al. (2002). (/) from Evans (2000). . 



dieted Cepheid masses larger than those derived from pulsation 
models (Neilson et al. 2011; Keller 2008; Bono et al. 2006, 
1999). To investigate the origin of this discrepancy, the combina- 
tion of spectroscopic and interferometric measurements will al- 
low us to derive orbital elements and dynamical masses. Recent 
analyses for two eclipsing binary Cepheids in the LMC were car- 
ried out by Pietrzynski et al. (2010, 201 1) which yielded masses 
with an accuracy < 1.5%. These dynamical masses are in agree- 
ment with the ones calculated from stellar pulsation models, 
suggesting that the pulsational theory provides the true current 
masses (see also Prada Moroni et al. 2012). However, this first 
conclusion needs additional mass estimates to better constrain 
the two models, particularly for wide binaries for which no sig- 
nificant physical interaction between the stars is expected. 

A number of Galactic Cepheids are known to have compan- 
ions closer than 30" (e.g. Evans et al. 2008; Evans 2011), but 
most of them are located too close to the Cepheid (~ 1-20 mas) 
to be observed with a 10-meter class telescopes. The already ex- 
isting orbit measurements were estimated only from IUE spec- 
trum or from the radial velocity variations. The only actual way 
to spatially resolve such systems is to use long-baseline inter- 
ferometry. We therefore have started a long-term interferomet- 
ric observing program that aims at studying a sample of seven 
northern and southern binary Cepheids. The first goal is to de- 
termine the angular separation and the apparent brightness ratio 
from the interferometric visibility and closure phase measure- 
ments. Our long-term objective, which needs a good sampling 
of the orbital period to get a reliable fit (several years), is to 
determine the full set of orbital elements, absolute masses and 
geometric distances. 

The binary (or multiple) Cepheid systems were selected ac- 
cording to two main criteria. Firstly, the angular separations 
(mainly estimated from spectroscopy so far) have to be resolv- 
able by the existing long-baseline interferometers. Secondly, the 
contrast between the Cepheid and its companion should not be 
too large, to ensure the flux contribution of the companion is 
detectable in the data. We therefore selected systems with an an- 
gular separation > 0.5 mas, and a dynamic range > 1 : 300. 

In this paper, we present the first results for the Cepheid 
VI 334 Cyg, observed with the MIRC combiner in the H band at 
the CHARA Array. The paper is organized as follows. We first 
introduce the observed Cepheid, with some knowledge about the 
primary and its companion. In Sect. 3, the instrument configura- 
tions and the data reduction are detailed. The data analysis, in- 
cluding the interferometric models used and the fitting steps, is 
discussed in Sect. 4. We then combined our interferometric re- 
sults with spectroscopic data in Sect. 5 to derive the parameters 
of the V1334 Cyg system. 



2. V1334 Cyg 

V1334 Cyg (HD 203156, HR 8157) was the first Cepheid ob- 
served for our program. We present in Table 1 the known param- 
eters of this system from the literature. The information about 
the Cepheid companion were derived from spectroscopic obser- 
vations (Evans 1995, 2000). 

This short-period Cepheid is an interesting system because 
it has been studied for many years but the binary or triple na- 
ture of the system is still debated. It has been suspected to be 
a member of a visual binary system for decades (see e.g. Millis 
1969; Abt & Levy 1970), with a separation between 0.1-0.2", 
but it has not been spatially resolved so far. Early radial ve- 
locity measurements also showed strong evidence of a spec- 
troscopic binary, but the orbital period could not be fully con- 
strained. Abt & Levy (1970) found a period of ~ 30 years for 
the visual component, but they also noticed a variation in the 
centre of mass velocity (v r ) of a shorter time-scale. The same 
variation was also observed by Szabados (1991), and is likely 
linked to a second closer companion. From the International 
Ultraviolet Exporer (IUE) low -resolution spectra, Evans (1995) 
detected the hottest star in the system, and derived the spectral 
type to be a B7.0V star. In a subsequent work, Evans (2000) 
concluded that the hottest star is the visual companion. From 
30 years of radial velocity measurements, the same author solved 
for the orbital parameters of the companion and found a period 
of 1937 days (see Table 1), that is significantly shorter than the 
value derived by Abt & Levy (1970), strongly suggesting the 
presence of a third component. They also derived a projected 
semi-major axis for the orbit of the Cepheid around the center 
of mass ai sin/ = 2.46 AU (see Table 1). Kiss & Vinko (2000) 
also detected the change in v y , and hypothesized a yellow-bright 
close companion. 

Many authors attempted to resolve the wide component, but 
we still do not have a firm conclusion because of its intermit- 
tent detection. For instance, Evans et al. (2006) set an upper 
limit of ~ 20 mas using the Hubble Space Telescope in the far- 
ultraviolet (although depending on the brightness of the stars). 
Scardia et al. (2008) mentioned the detection of the compan- 
ion from speckle interferometry in the V band, and measured a 
separation of 160 mas, while most of previous speckle observa- 
tions failed to resolve the system (from 1976 to 2005, see e.g. 
McAlister 1978; Hartkopf & McAlister 1984; Evans et al. 2006, 
see the last reference for a more detailled discussion), setting a 
separation < 35 mas. As argued by Evans et al. (2006), these in- 
termittent detections could be linked to a peculiar orbit, making 
the visual companion undetectable at some orbital phase, or it 
could even not exist. 
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Fig. 1. (u, v) plane coverage of our observations of VI 334 Cyg. 



3. Observations and data reduction 

Our observations were performed using the Michigan InfraRed 
Combiner (MIRC), installed at the CHARA array (ten 
Brummelaar et al. 2005) located on Mount Wilson, California. 
The array consists of six 1 m aperture telescopes with an Y- 
shaped configuration (two telescopes on each branch), oriented 
to the east (El, E2), west (W1,W2) and south (SI, S2), and so 
offering a good coverage of the (u, v) plane. The baselines range 
from 34 m to 331m, providing an angular resolution down to 
0.5 mas in H. 

MIRC (Monnier et al. 2004, 2010) combines the light com- 
ing from all six telescopes in the K or H bands, with three spec- 
tral resolutions (R - 42, 150 and 400). The recombination of 
six telescopes gives simultaneously 15 fringe visibilities and 20 
closure phase measurements, that are our primary observables. 

Our observations were carried out on UT 2012 July 27 and 
October 1 using an H band filter (Aq - 1 .65 fim) and either four, 
five or six telescopes. We used the low spectral resolution mode, 
where a prism splits the light on the detector into 8 narrow spec- 
tral channels. The (u, v) plane coverage for these nights is pre- 
sented in Fig. 1. We followed a standard observing procedure, i.e 
we monitored the interferometric transfer function by observing 
a calibrator before and after our Cepheids. The calibrators were 
selected using the SearchCal 1 software (Bonneau et al. 2006, 
201 1) provided by the JMMC. The journal of the observations is 
presented in Table 2, and the corresponding calibrators are listed 
in Table 3. 

The data were reduced with the standard MIRC pipeline 
(Monnier et al. 2007). The main procedure is to compute squared 
visibilities and triple products for each baseline and spectral 
channel, and to correct for photon and readout noises. A recent 
upgrade includes a simultaneous measurement of the photomet- 
ric channels with the fringes, enhancing the accuracy of MIRC 
down to 3 % in visibilities (Che et al. 2010, 2012). For the July 
observations, we used a coherent integration of 75 ms to improve 
the signal-to-noise ratio of the closure phase, at some expense of 



Table 2. Journal of the observations. 



UT 


Star 


Configuration 


2012 Jul. 27 09:52 


HD 200577 


S1-S2-E1-E2-W2 


10:23 


V1334 Cyg 


S1-S2-E1-E2-W2 


11:08 


HD 214200 


S1-S2-E1-E2-W2 


11:44 


V1334 Cyg 


S1-S2-E1-E2 


2012 Oct. 01 02:50 


HD 185395 


S1-S2-E1-E2-W1-W2 


03:40 


V1334 Cyg 


S1-S2-E1-E2-W2 


05:01 


HD 199956 


S1-S2-E1-E2-W1-W2 


05:44 


V1334 Cyg 


S1-S2-E1-E2-W1-W2 


06:28 


HD 218470 


S1-S2-E1-E2-W1-W2 


07:11 


V1334 Cyg 


S1-S2-E1-E2 


07:48 


HD 207978 


S1-S2-E1-E2 



Table 3. Calibrators used for our observations. 



Calibrator 


my 


w H 


Sp. Type 




r 


(HD) 








(mas) 


O 


200577 


6.1 


4.1 


G8ni 


0.758 ± 0.052 


3.2 


214200" 


6.1 


4.2 


K0 


0.790 ± 0.050 


15.5 


185395° 


4.5 


3.7 


F4V 


0.750 ± 0.060 


18.2 


199956 


6.6 


4.5 


K0 


0.603 ± 0.043 


7.1 


218470 


5.7 


4.7 


F5V 


0.477 ± 0.033 


15.5 


207978 


5.5 


4.4 


F6IV 


0.571 ±0.040 


21.6 



Available at http://www.jmmc.fr/searchcal. 



Notes, my, mu. magnitudes in V and H bands. 6>ud: uniform disk angu- 
lar diameter in H band, y: angular distance to the Cepheid. 
(n) Redetermined from internal MIRC calibration with HD 200577. 



the visibility calibration, but for October observations, we used a 
more standard 17 ms integration time. We applied the same cal- 
ibration error model as described in Monnier et al. (2012). We 
then did an incoherent average of 15 min for the final data. We 
note that this average prevent us from detecting a periodic sig- 
nal from a potential third component, and the detection of any 
incoherent light would then be limited by the uncertainty in the 
visibilities. 



3 



A. Gallenne et al.: Multiplicity of Galactic Cepheids from long-baseline interferometry 



4. Model fitting 

To model the squared visibilities, triple amplitude and clo- 
sure phase signals, we used the LITpro 2 model fitting software 
(Tallon-Bosc et al. 2008), based on the Levenberg-Marquardt al- 
gorithm. It provides a set of elementary models that can be com- 
bined all together. The software also contains a tool allowing 
the search for the global minimum to solve for the problem of 
multiple x 1 minima. 

The two epoch observations (July and October) were re- 
duced separately in order to detect the changing position of the 
companion. 

4.1. The models 

As a first step, the primary component (the Cepheid) was mod- 
eled with a uniform disk (UD) angular diameter. The complex 
visibility model is: 

V*(m,v) = , 



Table 4. Summary of the parameters estimated from the model 
fit. 



with J\{x) the first-order Bessel function, x = 7r6> UD Vm 2 + v 2 /A, 
(u, v) the spatial frequencies, 6>ud the UD angular diameter, and 
A the wavelength. 

The choice of a UD diameter instead of a limb-darkened 
(LD) disk for the fitting procedure is justified because the angu- 
lar diameters of the Cepheids are small compared to the angular 
resolution of the interferometer, and the limb darkening effects 
are therefore undetectable. The conversion from UD to LD an- 
gular diameter was done by using a linear-law parametrization 
h(p) — 1 _ M i(l _ A*X with the LD coefficient u A from Claret & 
Bloemen (2011). The conversion is then given by the approxi- 
mation (Hanbury Brown et al. 1974): 



1-7^/15 



It is worth mentioning that the uncertainty of the limb- 
darkened coefficient has a small impact on the angular diame- 
ter conversion. A variation of 20 % of ux gives a LD diameter 
difference of less than 0.5 %. 

As a second step, we fitted to the data a model with one com- 
panion, assumed to be unresolved by the array. The correspond- 
ing model for this binary system is: 

y*( M ,v) + /y(M,v) 



V(u,v) 



where / is the flux ratio between the companion and the 
Cepheid, and V(u, v) is the complex visibility model of an un- 
resolved source: 

V(u, v) = exp(-2in(uAa + vA6)/A), 

with (Aa, A6) the relative position of the companion w.r.t the 
Cepheid. 

The closure phase is then estimated from the modulus and ar- 
gument of the bispectrum, B123, for each closed baseline triangle 
and spectral channel: 

|Bl2 3 | = \V(m,Vl)V(u 2 ,V2)V*(U3,V 3 )\ 

arg(B I2 3) = 0123 = arg(y(Mi,vi)y(M 2 ,V2)V*(«3, v 3 )) 

We assumed that the variation of the angular diameter be- 
tween different acquisitions in the same night is negligible com- 
pared to our level of accuracy. 





2012-07-27 


2012-10-01 


Single star model 






e UD (mas) 


0.565 ± 0.052 


0.487 ± 0.045 


9 LD (mas) 


0.575 ± 0.052 


f\ A C\/~ f\ f\ A C 

0.496 ± 0.045 




1.63 


2.08 


Binary model 
0ud (mas) 


0.494 ± 0.053 


0.436 ± 0.045 


9 LD (mas) 


0.503 ± 0.053 


0.444 ± 0.045 


/(%) 


3.15 ± 0.15 


3.08 ± 0.09 


Aa (mas) 


-1.153 ±0.030 


-0.113 ±0.014 


A6 (mas) 

2 

XI- 


-8.836 ± 0.017 


-8.359 ± 0.009 


0.34 


1.24 



LITpro software available at http://www.jmmc.fr/litpro. 



Notes. 6 VD and 9 LD are the uniform and limb-darkened disk angular 
diameter, respectively. /, x and y correspond to the flux ratio and posi- 
tion of the companion. x 1 is the reduced x 1 of the corresponding fitted 
model. 



4.2. Fitting steps 

We performed a least-squares model fit simultaneously with the 
squared visibility, triple amplitude and closure phase measure- 
ments. Our search strategy was the following. We first proceeded 
to a grid search in the^- 2 space to determine the approximate po- 
sition and brightness ratio of the companion. A first grid search 
between +20 mas with a 0.2 mas spacing was performed with 
various flux ratio (from < / < 0.2 with 0.005 steps). Then a 
second grid of +1 mas with a 0.01 mas spacing around the most 
likely position and brightness ratio was used to refine the posi- 
tion. Finally all parameters were fitted using the refined values. 
The first guess for the angular diameter was taken from Table 1 . 

4.3. Results 

The closure phase signal and squared visibilities, presented in 
Figs. 2, 3 and 4, clearly show a departure from a single star with 
a symmetric brightness distribution. These variations reveal the 
presence of at least one companion. 

We chose as first guess 6*ud = 0.547 mas, and performed the 
fitting procedure with a grid search as explained in Sect. 4.2. The 
probability maps for our two epochs of observations are shown 
in Fig. 5, and the fitted parameters are reported in Table 4. The 
model is also represented graphically with black dots in Figs. 2, 
3 and red dots in Fig. 4. We notice a good agreement between 
the model and the data. The fitted model give the most probable 
location of this companion at an angular separation p = 8.91 mas 
and a position angle PA = -172.6° for the first epoch, and p = 
8.36 mas and PA = -179.2° for the second. The measured flux 
ratio is also particularly consistent between the two epochs. 

We estimated the limb-darkened angular diameter to be 
#ld = 0.503 ± 0.053 mas for the first epoch, using the limb- 
darkening coefficient ux = 0.2423 (Claret & Bloemen 2011), 
chosen from the stellar parameters: T e ff = 6250 K, logg = 2.0, 
[Fe/H] = 0.0 and v t = 4 m s" 1 (Luck et al. 2008, at pul- 
sation phase = 0.17). For the second epoch, the photo- 
spheric parameters did not change significantly (because it is 
low-amplitude pulsation Cepheid), leading to a same u A to find 
#ld = 0.444 ± 0.045 mas (at pulsation phase 4> = 0.91). These 
measurements are in agreement with the averaged angular diam- 
eters predicted from a surface brightness relation (~ 0.545 mas, 
Kervella et al. 2004b; Moskalik & Gorynya 2005), taking into 
account the small amplitude variation. As no IR photometric 
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A (^m) A (^m) A (^m) A (^m) A (/im) 

Fig. 2. Closure phase signal for the first epoch. The color-coded squares are the data, while the black dots represent the binary 
model. 




A (^m) A (^m) A (^m) A (^m) A (/im) 



Fig. 3. Closure phase signal for the second epoch. The color-coded squares are the data, while the black dots represent the binary 
model. 

curves are available, it was not possible to compare our value 5. Combination of spectroscopy and interferometry 

with an angular diameter derived from a SB technique for the 

actual phase of our measurements . 5.1. Preliminary orbit 

Because of its intermittent detections, the existence of the wide 
component is rather uncertain, and might not exist. We claim that 
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Fig. 4. Squared visibility measurements for the two epochs. The data are represented in blue, while the red dots are the fitted binary 
model. 
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Fig. 5. Probability maps for the two observed epochs, showing the maximum obtained at each point of the search region with a 3 % 
flux ratio. 



the companion detected with MIRC and IUE (Evans 1995) can 
be identified as the spectroscopic one. Therefore, our measured 
separations should be compatible with the spectroscopic orbit, 
and it should be possible to derive all the orbital parameters. 

We therefore combined our astrometric measurements with 
the single-line radial velocity data gathered by Evans (2000), in 
order to derive the complete orbital elements for VI 334 Cyg Ab. 
We stress that this analysis is preliminary as we have only two 
astrometric measurements so far. We computed the orbital so- 
lutions through the formalism developed by Wright & Howard 
(2009), slightly modified to include the pulsation velocity of the 



form: 



Vpui s (0 = J]Ai cos(2/7r^(f) + BO 



where A{ and fi; are the fitted Fourier coefficients, and <p denotes 
the pulsation phase. 

We solved for all the orbital elements 
(Porb, T v ,e, K, v y , a>,a, i, Q) and pulsation parameters 
{Ppuk,A\,A2,Bi,B2). The reference epoch of maximum 
light for the Cepheid, Tq, was held fixed to the value given by 
Samus et al. (2009). The initial values for the orbit and the 
pulsation were chosen from Evans (2000) and Samus et al. 
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Fig. 6. Left : fitted (solid line) and measured (blue dots) orbital velocity. Middle: fitted (solid line) and measured (blue dots) pulsation 
velocity. Right: orbit of V1334 Cyg Ab. The data points are the MIRC results from Table 4. 



(2009). The final elements derived from our combined fit are 
listed in Table 5. The fitted parameters are in good agreement 
with those from Evans (2000). The quoted uncertainties for the 
elements derived from spectroscopy (P 0{ \„T v ,e, K, v y ,a>,) were 
estimated using the bootstrapping technique (with replacement 
and 500 bootstrap samples). For the remaining elements derived 
from interferometry (a, i and Q), we refitted the orbits 500 
times, each time adding Gaussian noise to each astrometric 
point according to their uncertainties. The standard deviation 
from these trials is then used as the uncertainty. Fig. 6 (left and 
middle) shows the orbital and pulsation velocities disentangled 
from the radial velocity measurements. The solid black lines 
denote our fitted curves. The final best fit orbit of VI 334 Cyg 
Ab is also plotted in Fig. 6 (right) with our MIRC measurements 
marked by the blue dot symbols. 



5.2. Apparent magnitude, spectral type and mass of the 
companion 

Combining the //-band magnitude toh = 4.66 + 0.04 given by 
the 2MASS catalog (Cutri et al. 2003) with our averaged mea- 
sured flux ratio / = 3.10 + 0.08 %, we derived for the com- 
panion a magnitude mn(comp) = 8.47 + 0.15, and mn(cep) 
4.70 + 0.15 mag. As no //-band light curve is available to es- 
timate the Cepheid magnitude at our pulsation phase, an addi- 
tional uncertainty of 3 % was quadratically added to take into 
account the phase mismatch. The choice of these 3 % is based 
on the amplitude variation of the light curve in V (Klagyivik & 
Szabados 2009, that is surely to be lower in H). 

The absolute magnitude, Mh, can be estimated knowing the 
distance to the system. However, there is no accurate determina- 
tion of the distance for this Cepheid. The Hipparcos data give a 
distance d = 662 + 162 pc (van Leeuwen 2007). The use of a K- 
band P-L relation for first overtone (FO) pulsators (Bono et al. 
2002, non-canonical model) gives d = 683 ± 17 pc, while con- 
verting the overtone period to the fundamental one with the pe- 
riod ratio data from Alcock et al. (1995) and using a Zf-band P-L 
relation for fundamental (F) mode pulsator (Storm et al. 2011) 
gives d = 639 + 17 pc. As there is no optimum value, we plotted 
in Fig. 7 (top panel) the spectral type vs. the distance. The pre- 
vious cited distance ranges set a spectral type for the companion 
between a B8.0V and B4.0V star. The extinction was assumed 
negligible at our observing wavelength (Ah = 0.023 mag, es- 



Table 5. Fitted orbital and pulsation parameters of VI 334 Cyg 
Ab. 





Spectroscopy only 


This work 




(Evans 2000) 




Orbit 






Porb (days) 


1937.5 ±2.1 


1938.6 ± 1.2 


T p (HJD) 


2 443 607 ± 14 


2 443 616.1 ±7.3 


e 


0.197 ±0.009 


0.190 ±0.013 


Ki (km s" 1 ) 


14.1 ±0.1 


13.86 ±0.17 


v y (km s~') 


-1.8 ±0.1 


-1.9 ± 0.1 




226.3 ± 2.9 


228.7 ± 1.6 


no 




206.3 ± 9.4 


a (mas) 




8.54 ±0.51 


((°) 




124.7 ± 1.8 


m H 




8.47 ±0.15 


Pulsation 






-Ppuis (days) 


3.33251 ±0.00001 


3.33250 ± 0.00002 


T a (HJD) 


2 440 124.5330 


2 440 124.5330 


A, 




4.35 ±0.15 


A 2 




1.81 ±0.11 


Si 




0.08 ± 0.06 


B 2 




2.72 ± 1.30 



Notes. P 0I t,: orbital period. T p : time passage through periastron. 
centricity. K: radial velocity semi-amplitude of the primary. v y 



>: ec- 
sys- 



temic velocity, oj: argument of periastron. Q: position angle of the as- 
cending node, a: semi-major axis, i: orbital inclination. m H : apparent 
magnitude in H. P pu \ s : pulsation period. T : reference epoch of maxi- 
mum light. A,, B\\ Fourier parameters. 
<n) From Samus et al. (2009), and held fixed when fitting 



timated using the total-to-selective absorption ratios Ry = 3.1, 
Ra = An/E(B - V) = /?v/6.82 from Fouque et al. 2003, and 
the average color excess E(B - V) = 0.05 from Evans 1995 and 
Kovtyukh et al. 2008) 

We can also estimate the mass ratio, q, for a range of dis- 
tances by combining the Kepler's third law with the spectro- 
scopic mass function, 



M T - Mi + M 2 - 



a 3 d 3 

^orb 



(M 2 sin 3 
(Mi + M 2 ) 2 



3.784 x 10~ 5 K 3 P, 



1 ^orb 



2n3/2 



(1) 

(2) 
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450 500 550 600 650 700 750 800 850 



d (pc) 

Fig. 7. Absolute magnitude, spectral type of the companion (top 
panel), and mass of each component (bottom panel) as function 
of the distance. The shaded gray areas denote the distance range 
of Hipparcos, F mode P-L relation (left), and FO mode P-L re- 
lation (right). The vertical dotted line represents the minimum 
distance where Mi = Mi- 



which yields: 



M2 



a d sin i 



0.03357tfiP orb Vl -e 1 



with a in arcseconds, d in pc, P OI \, in yr, and K t in kms l . 
Combining Eq. 1 and 2 gives the mass of each component: 



M 2 = 0.03357 



K,a 2 d 2 Vl -e 2 



P 0lb sin i 



M x = M T - M 2 . 



(3) 
(4) 



We plotted Eq. 3 and 4 vs. the distance in Fig. 7 (bot- 
tom panel). The shaded gray areas denote the distance range of 
Hipparcos, F mode P-L relation, and FO mode P-L relation. As 
we expect a lower mass companion (Evans 1995), we can set a 
lower limit to the distance, q < 1 implies d > 691 pc, plotted as a 
dotted vertical line in Fig. 7. This limit is consistent with the dis- 
tance derived from the FO P-L relation (d = 683 pc, Bono et al. 
2002), and would mean a mass ratio equal to 1 for this system. 
With this limit, Eq. 3 yields a minimum mass for the companion 
(and so the Cepheid) — 3.6 M B , and a spectral type earlier 
than a B5.5V star. 



6. Conclusion 

We presented new multiple telescope interferometric observa- 
tions of the classical Cepheid V1334 Cyg in the H band. For 
the first time, we were able to spatially resolve the companion at 
two epochs. We derived the limb-darkened angular diameter for 
the Cepheid at their corresponding pulsation phase, the relative 
positions, and the flux ratio of the companion. We combined our 
accurate astrometric measurements with existing spectroscopic 
data and derived preliminary orbital solutions for the VI 334 Cyg 
system. We also determined a minimal distance to the system to 
be d ~ 691 pc, which yields to a minimum mass for the compo- 
nent M2 ~ 3.6 M . We also found that its spectral type is earlier 
than a B5.5V star. 



Our work, using multi-telescope recombination, provided 
unique and useful informations, both on the Cepheid and its 
companion. These innovative results show the capabilities of 
long-baseline interferometry to study Cepheids in binary or mul- 
tiple systems. This is particularly important as most of the com- 
panions are located too close to the star to be spatially resolved 
by a single telescope. This technique provides new observables 
that can be efficiently combined with spectroscopic results to 
provide innovative constraints on the system properties. We plan 
to continue our observing program in the next years to increase 
the sample of binary Cepheids with well determined orbital ele- 
ments. 
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